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ABSTRACT 

Context. The stellar system S Vel contains the brightest eclipsing binary in the southern sky {d Vel A), and a nearby third star located 
a 0.6" away ((5 Vel B). The proximity of (5 Vel B (usable as a reference) makes it a particularly well-suited target to detect the 
astrometric displacement of the center of light of the eclipsing pair 

Aims. We obtained NACO astrometric observations with two goals: 1) to confirm the orientation of the orbital plane of the eclipsing 
pair on the sky determined by interferometry (Paper III), and 2) to demonstrate the capabilities of narrow-angle adaptive optics 
astrometry on a simple system with predictable astrometric properties. 

Methods. We measured the angular separation vector between the eclipsing binary 6 Vel A and the visual companion S Vel B from 
narrow-band images at 2.17 //m obtained with the VLT/NACO adaptive optics system. Based on these observations and our previous 
determination of the orbital parameters of the wide binary 5 Vel A-B, we derived the apparent displacement of the center of light of 
the eclipsing pair at 1 1 epochs over its orbital cycle. 

Results. We detect the astrometric wobble of the center of light of the 6 Vel A pair relative to B with a typical measurement precision 
of X 50 fias per epoch, for a total amplitude of the measured displacement of a; 2 mas. 

Conclusions. The detected wobble is in relatively good agreement with the model we presented in Paper III, and confirms the 
orientation of the Aab orbital plane on the sky. The residual dispersion compared to our model is 110/jas rms, which we tentatively 
attribute to photometric variability of the fast rotating A-type components Aa and/or Ab in the Bry line. Based on these results, 
we conclude that in favorable conditions (bright source with only two resolved components, small angular separation), narrow-angle 
astrometry with adaptive optics on an 8-meter class telescope can reach an accuracy of 50 to lOOyuas. 

Key words. Stars: individual: (HD 74956, 6 Vel); Stars: binaries: eclipsing; Methods: observational; Techniques: high angular 
resolution; Astrometry 



1. Introduction 

The bright southern star 6 Vel (HD 74956, HIP 41913, GJ 321.3, 
GJ 9278) is a multiple system comprising at least three stars. Its 
brightest component, 6 Vel A, was identified only in 2000 as one 
of the brightest eclipsing system in the sky ( [Otero et al.|2000j l. 
We consider here the two main components of the 6 Vel system 
(Aa-i-Ab and B) as a test case to validate the potential of high- 
precision narrow-angle astrometry from the ground using adap- 
tive optics. Exciting results were recently obtained using this 
observing technique, e.g., by Clarkson et al. (2012 1, |Hufimann^ 
^TaLl ( |2OT2l i, and [Kolher, Ratzka & Leinertr ( |2012| i. These au- 
thors typically measured differential astrometry within a field of 
view of a few arcseconds to an accuracy of a few 100/ias (see 
e.g. Table 4 in |Clarkson et aL|2012[ ). We propose here to push 
this observing technique to its precision limits on 5 Vel, which 
appears to be a particularly favorable target (bright and geomet- 
rically simple), in order to evaluate its potential. 



In previous works on the 6 Vel system, Kervella et al. ( 2009 
Paper I) searched for infrared circumstellar excess in the irmer 
6 Vel A-B system, whose present angular separation is approx- 
imately 0.6". Part of the NACO data presented here was also 
used in this first article, although not for high-precision astrom- 
etry. IPribulla et aT] (jSUTTI Paper II) and [Merand et aT] ( [2011] 
Paper III) obtained the orbital elements and fundamental param- 
eters (including the parallax of the system) of the three compo- 
nents of the 6 Vel system. In the present study, we take advan- 
tage of the well-constrained and accurate model of the system 
obtained in these previous works to predict the astrometric wob- 
ble of the center of light (hereafter CL) of the eclipsing system, 
and compare it with the measured NACO astrometry. We detail 
our observations and data analysis procedure in Sect. |2] and de- 
scribe the method we used to derive the wobble of the CL of the 
eclipsing component in Sect. [3] We discuss in Sect. [4] the dif- 
ferences between the astrometric wobble predicted by the model 
and the observed displacement. 



Send offprint requests to: P. Kervella 

* Based on observations made with ESO telescopes at the La Silla 
Paranal Observatory, under ESO programs 08 l.D-0 109(A), 081.D- 
0109(B), 382.D-0065(A), and 282.D-5006(A). 
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2. Observations and data reduction 

2.1. Observing log and raw data processing 

We observed 6 Vel at 1 1 epochs in April-May 2008 and January 
2009 using the Nasmyth Adaptive Optics System (NAOS, Rous- 



|set et al. |[2003) of the Very Large Telescope (VLT), coupled to 
the CONICA infrared camera ( |Lenzen et al.|1998| l, which is ab- 
breviated as NACO. These observations were obtained primarily 
to provide high-precision differential astrometry of the eclipsing 
system 5 Vel A relative to B, following the promising demon- 
stration bylSeifah rt, Roll & Neuha user (2008). We selected the 
smallest available pixel scale of 13 mas/pixel, giving a field of 
view of 13.6"xl3.6". Because of the brightness of 5 Vel, we em- 
ployed a narrow -band filte:[^at a wavelength of 2.166 + 0.023 //m 
(hereafter abbreviated as 2.17yum) together with a neutral den- 
sity filter (labeled "ND2_short", that has a transmission of about 
1.5%). The 2.17/im filter wavelength corresponds to the Brack- 
ett 7 line of hydrogen. A detailed list of the NACO exposures ob- 
tained in 2008 is presented in Table 2 of Paper I. At each epoch 
50 exposures were obtained, always centering component A on 
the same pixel of the CONICA detector array to mitigate the 
possible effects of field distorsion. 

The raw images were processed in a standard way: dark 
subtraction, flat fielding, and cosmetics removal. We consid- 
ered three approaches for the flat-field and dark-frame calibra- 
tion: using specific dark and flat-field frames obtained sepa- 
rately for each night, using average flat/dark frames, and using 
median flat/dark frames (the last two computed over our 2008 
observation epochs). For these different methods, the star posi- 
tions extracted from the resulting data cubes are identical within 
4 X 10" ■'pixels, or 50juas. This figure is comparable to our sta- 
tistical measurement uncertainties, and we conclude that we do 
not detect a significant astrometric variability induced by the flat 
field or dark. We therefore chose to use the median flat and dark 
frames over all observations. This remarkable stability is proba- 
bly because NACO was kept cold and was not dismounted from 
the telescope over the period of our 2008 observations. Along 
this line, Neuhaus er et al.| ([2008 ) demonstrated that the plate 
scale of NACO is extremely stable (to < 5 x 10""*, in relative 
terms) over three years. The cosmetic corrections we applied 
concerned only pixels located outside of our astrometric win- 
dows (Sect. 2.2 1, as the areas used for each star position were 
specifically chosen to present no dead or hot pixels. They there- 
fore had no influence on the astrometric measurements. 



2.2. Star separation measurement 

To measure the separation of the two star images on the NACO 
detector, we treated each image in our data cubes separately. We 
extracted two sub-images of 21 x 21 pixels centered on the posi- 
tions of A and B (as shown in Fig.[T]). We then used a classical 
minimization algorithm to shift and scale the sub-image of 
(5 Vel A over the position of the fainter component B. Our choice 
of relatively small sub-windows prevents the contamination by 
the background noise that may present inhomogeneous patterns. 
We checked that the fit is not sensitive to the chosen size of the 
sub-windows over a large range of values. To achieve this fit, 
we computed a sub-pixel interpolation of the shifted image of A, 
using the Yoricl^language, and more specifically using an inter- 
polation routine in the Fourier spacej^ We derived four parame- 

' http://www.eso.org/sci/facilities/paranal/instruments/naco/ 
^ http ://y orick. sourcef orge . net/ 

f ft_f ine_shift function of the fft_utils . i package 



Fig. 1. Average image of 6 Vel obtained on 2008 Apr 01, together with 
the 21x21 pix astrometric windows that we used for the differential 
astrometry (white squares). 



ters: the relative shifts dx and dy in pixels between the two star 
images, the flux ratio p = fAlfa, and the uniform background 
level C. In order to estimate their associated error bars, we used 
the bootstrapping technique described in Appendix B of Ke rvella| 
et al. ( 2004j ). We validated the Fourier interpolation method by 
comparing its results with a simple Gaussian fit of the two PSF 
cores that give the same relative positions (within 150/ias) al- 
though with larger dispersion due to the mismatch of the slightly 
seeing-distorted PSF and the Gaussian function. Thanks to their 
angular proximity, 5 Vel A and B share the same residual opti- 
cal aberrations from the adaptive optics system. Because it takes 
advantage of the similarity of their images, our fitting approach 
based on Fourier interpolation is particularly efficient. 

For the conversion of the star separations measured in pix- 
els to angular values, we adopted the pixel scale of 13.26 + 
0.03 mas/pixel (Masciadri et al. 2003). This value is in good 



agreement with the pixel scale obtained by Neuhauser et al. 
( |2008) 13.24 ± 0.06 mas/pix). Our differential astrometric mea- 
surements are affected only by the +0.23% pixel scale uncer- 
tainty over the amplitude of the astrometric wobble of the CL 
of Aa-Ab and not over the full angular separation of A and B. 
As shown in Sect. [3] the astrometric wobble amplitude is less 
than 2 mas, and the corresponding ^ 4/ias systematic error is 
therefore negligible compared to our statistical accuracy. 



2.3. Differential atmosplneric refraction 

The angular separation of S Vel A and B is only « 0.6" for the 
epoch of our observations. Although this angle is quite small, 
the effect of differential atmospheric refraction is not negligible 
considering that the precision of our measurements (s; 50 yuas per 
epoch in right ascension and declination) corresponds to only 
~ 1/12 000* of the separation of the two stars. Figure |2] shows 
schematically the effect of atmospheric refraction on the appar- 
ent positions of the two stars. In this diagram the sky positions 
of the two stars in the absence of atmosphere are labeled with 
the "vac " subscript, while the apparent positions measured from 
the ground are labeled "atm"; tp is the parallactic angle (positive 
after local celestial meridian crossing) at the time of the obser- 
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Fig. 2. Geometrical effect of the atmospheric refraction on the differ- 
ential astrometry of a binary star (see Sect.|2.3|for details). 



Table 2. Photometric flux ratio /(A)//(B) 
NACO 2. 17yum narrow-band filter. 



between A and B in the 



UT date 


MJD (f> 


/(A)//(B) 




-54000 





2008 Apr 01 


557 


0224 


0.9294 


10.547 


+ 


0.009 


2008 Apr 04 


560 


9976 


0.0175 


10.413 


+ 


0.012 


2008 Apr 06 


562 


0121 


0.0400 


10.426 


+ 


0.008 


2008 Apr 07 


563 


0048 


0.0619 


10.332 


+ 


0.007 


2008 Apr 20 


576 


9715 


0.3713 


10.469 


+ 


0.008 


2008 Apr 23 


579 


0231 


0.4167 


10.147 


+ 


0.012 


2008 Apr 24 


580 


9917 


0.4603 


10.589 


+ 


0.011 


2008 May 05 


591 


9748 


0.7036 


10.295 


+ 


0.011 


2008 May 07 


593 


9732 


0.7478 


10.269 


+ 


0.011 


2008 May 18 


604 


0442 


0.9709 


6.993 


+ 


0.009 


2009 Jan 07 


838 


1347 


0.1556 


10.385 


+ 


0.022 



Notes. The 2008 May 18 epoch corresponds to the primary eclipse. 



vation; Vatm and Vvac are the vector separations of A and B with 
and without atmosphere, respectively; and are the apparent 
vector displacements of the position of the stars due to atmo- 
spheric refraction; and r is the resulting differential astrometry 
correction vector. 

The quantity of interest for our purpose is the vector separa- 
tion of the two stars Vvac ( Aff vac > ^^vac)- From a simple vectorial 
closure, the transformation between the measured differential as- 
trometry vector on the sky Vatm (Aff vac, A^atm) is 

Vvac = Vatm + r, (1) 

where r = eA ~ fB- We computed the values of IcaI and |eB| us- 
ing the slaRefro routine from the Starlink librar}|^and the 
apparent positions of the two stars. The r = (Aa, , A6r) vector 
is always pointing towards zenith. This routine realizes an inte- 
gration of the effect of the atmospheric crossing on the apparent 
direction of the incoming star light using a realistic atmospheric 
model ( |Liebe et al.[ 1993 ). As input parameters, we used the ac- 
tual conditions measured at Paranal for each observation epoch 
(airmass, pressure, temperature, humidity, etc.). Further details 
on this routine can be found in [Gubler & Tytler| ( |1998| l. The 
right ascension and declination corrections to add to the mea- 
sured NACO differential astrometry of B relative to A are 

Aur = \r\smip, A6r = \r\cos(p. (2) 

For an additional discussion of the effects of atmospheric refrac- 
tion on adaptive optics observations, the interested reader is re- 
ferred to .Roe, ( ,2002j and Helminiak, ( ,2008) . 

3. Apparent astrometric displacement of dVel A 

The measured relative separations between 5 Vel A and B, cor- 
rected for differential atmospheric refraction, are listed in Ta- 
ble [T] The photometric flux ratio of the two stars is presented 
in Table |2] The orbital phases of the measurements with re- 
spect to the eclipsing binary Aab were computed using the fol- 
lowing parameters from Paper II (ROCHE model): ro(HJD) - 
2452528.950 (epoch of periastron passage) and P = 45.15023 d. 
The observed displacement of the CL of Aa-Ab relative to B 
results from the vectorial sum of the orbital motion of B on its 
142-year orbit, and the displacement of the CL of the eclipsing 

'* http ://starlink.j ach.hawaii . edu/ 



binary system over its » 45-day orbital cycle. In order to mea- 
sure this second component only, the secular orbital displace- 
ment of 6 Vel B relative to A must be subtracted from the mea- 
sured differential astrometry. 

In Paper III |Merand et al.| ( |201 l| l determined the orbital el- 
ements of the 5 Vel A-B system. We used these parameters to 
estimate the relative positions of the center of mass of (5 Vel A 
and B at each observation epoch, and subtracted the resulting 
vectors from our NACO astrometric measurements. The evolu- 
tion of the measured astrometric position of 5 Vel B relative to 
A as a function of time (over more than a century) is presented 
in Table [3] The model positions for the different epochs are also 
listed, with the residuals compared to the observations. 

Because of uncertainties in the determination of the orbital 
parameters of the A-B system, we added a constant astromet- 
ric shift (Aa - -0.3 mas, A6 = mas) to our CL displacement 
model to match the measured CL positions. To determine this 
constant shift, we minimized the dispersion of the 2008 mea- 
surements with respect to the model. We needed this coiTec- 
tion, as our NACO measurements have a limited absolute accu- 
racy caused by the plate scale uncertainty (a: 2 x 10"^). For the 
0.6" separation of 6 Vel A and B, this systematic uncertainty is 
a! 1 mas and it does not average out with the number of mea- 
surements. We obtained a residual rms dispersion of the NACO 
measurement points of 1 10/ias compared to our model (exclud- 
ing the 2008 May 05, 2008 May 07, and 2009 Jan 07 epochs, see 
Sect.[4T]i. 

The measured astrometric wobble of the CL of the eclips- 
ing system is presented in Fig, p] and the corresponding numer- 
ical values are listed in Table I?] Another view of this apparent 
displacement, including a plot of the differential photometry of 
(5 Vel A-B, is shown in Fig.|5] In these two figures, the predic- 
tion from the physical model of the eclipsing pair presented in 
Paper III is also shown for reference. The derived orientation of 
the orbital plane of the eclipsing pair is consistent with the in- 
terferometric measurements reported in Paper III. The observed 
amplitude of the CL displacement is also in satisfactory agree- 
ment with the predicted value (a 2 mas). 

It is interesting to note that at the phases of the eclipses in the 
model curves, the CL of Aab moves over small arcs because of 
the slight inclination of the orbital plane on the line of sight and 
also because of the inhomogeneous light distribution on the ap- 
parent disks of the two stars (they are both relatively fast rotators 
and therefore show some gravity darkening). The amplitude of 
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Table 1. Differential astrometry of (5 Vel derived from the NACO images. 



UT date 


MJD 

-54000 


HJD 

-2454000 





6 


AM 


Aa,- 




A^, 


[as 


- Ca] ± CTstat 


[Sb 


- 6a]+ CTstat 


2008 Apr 01 


557.0224 


557.5245 


0.9294 


0.77 


1 


16 


-0.014 





113 


-430 


570 ± 0.039 


457 


819 + 0.043 


2008 Apr 04 


560.9976 


561.4996 


0.0175 


0.96 


1 


16 


-0.028 





098 


-430 


669 ± 0.043 


457 


070 + 0.055 


2008 Apr 06 


562.0121 


562.5141 


0.0400 


0.68 


1 


16 


-0.010 





115 


-430 


828 ± 0.028 


457 


046 + 0.037 


2008 Apr 07 


563.0048 


563.5067 


0.0619 


0.92 


1 


16 


-0.015 





111 


-430 


724 ± 0.028 


456 


975 + 0.038 


2008 Apr 20 


576.9715 


577.4732 


0.3713 


0.68 


1 


16 


-0.009 





113 


-430 


267 ± 0.035 


452 


865 + 0.039 


2008 Apr 23 


579.0231 


579.5247 


0.4167 


1.02 


1 


19 


0.038 





063 


-429 


964 ± 0.030 


451 


974 + 0.033 


2008 Apr 24 


580.9917 


581.4932 


0.4603 


1.30 


1 


16 


0.026 





095 


-429 


790 ± 0.042 


451 


090 + 0.056 


2008 May 05 


591.9748 


592.4761 


0.7036 


0.79 


1 


17 


0.035 





079 


-428 


888 ± 0.046 


448 


714 + 0.042 


2008 May 07 


593.9732 


594.4744 


0.7478 


0.86 


1 


18 


0.036 





073 


-429 


234 ± 0.036 


447 


959 + 0.033 


2008 May 18 


604.0442 


604.5450 


0.9709 


0.76 


1 


48 


-0.040 


-0 


013 


-429 


604 ± 0.055 


445 


438 + 0.051 


2009 Jan 07 


838.1347 


838.6357 


0.1556 


1.20 


1 


41 


-0.207 


0.086 


-424.775 ± 0.084 


381.374 + 0.116 



Notes, (p is the orbital phase of the system (relative to periastron passage), 6 the DIMM seeing in arcseconds (in the visible), AM the airmass, and 
Aa,. and Ad,- the corrections for the differential atmospheric refraction (see Sect. |2.3| for details). The listed differential astrometry [a^ - cr^] and 
[6b - Sa] include the atmospheric refraction correction. The angles are all expressed in milliarcseconds. 



Table 4. Astrometric displacement of the center-of-light (CL) of S Vel A. 



UT date 


MJD (P 


Model A-B 


Model CL 


Obsei-v. CL 


No- 




-54000 


vector 


displacement 


displacement 





2008 Apr 01 


557 


0224 


0.9294 


-430 


397,457 


959 


+0.195,-0 


390 


-0.173±o.039,-0.140±o.043 


10 


2008 Apr 04 


560 


9976 


0.0175 


-430 


355,456 


900 


-0.1 17, +0 


278 


-0.314^0.043, +0.170^0.055 


4 


2008 Apr 06 


562 


0121 


0.0400 


-430 


344,456 


630 


-0.191, +0 


434 


-0.484±o.o28,+0.416±o.o37 


8 


2008 Apr 07 


563 


0048 


0.0619 


-430 


334,456 


365 


-0.255, +0 


566 


-0.390±o.028,+0.610±o.038 


4 


2008 Apr 20 


576 


9715 


0.3713 


-430 


183,452 


639 


-0.117, +0 


201 


-0.084^0.035, +0.226±o.o39 


1 


2008 Apr 23 


579 


0231 


0.4167 


-430 


160,452 


091 


-0.051,-0 


009 


-0.196±o.030, -0.117^0.033 


8 


2008 Apr 24 


580 


9917 


0.4603 


-430 


139,451 


565 


+0.070, -0 


201 


+0.349±o.o42, -0.475±o.o56 


7 


2008 May 05 


591 


9748 


0.7036 


-430 


015,448 


627 


+0.449, -0 


984 


+ 1.127±o.046, +0.087±o.042 


28 


2008 May 07 


593 


9732 


0.7478 


-429 


992,448 


093 


+0.464,-1 


006 


+0.758±o.o36,-0.134±o.o33 


26 


2008 May 18 


604 


0442 


0.9709 


-429 


876,445 


399 


+0.052, -0 


082 


-0.272±o.o55,+0.039±o.05i 


5 



Notes. "Model A-B vector" is the predicted astrometric separation ([o-g - a^], [6b - 6a]) of the centers of light of A and B (with a shift of 
Aa = -0.3 mas, see Sect.[3]l. The shift of the center of light of Aab relative to its center of mass is listed in the column "Model CL displacement". 
The observed Aab center of light displacement ("Observ. CL displacemenf ) and the distance between the observed positions and the model 
positions in number of standard deviations (Ncr) are also listed. All angular values are expressed in milliarcseconds. 



this effect is too small, however, to be detectable in our NACO 
measurements. 

4. Discussion 

4. 1 . Inconsistent NACO measurements 

The overall agreement between our NACO astrometric measure- 
ments and the model of the eclipsing pair is satisfactory. How- 
ever, three of our measurement points show a significant depar- 
ture from the expected CL trajectory, at epochs 2008 May 05, 
2008 May 07, and 2009 Jan 07. 

For our 2009 epoch, the observed discrepancy is most prob- 
ably due to the combination of two causes. The first is the insta- 
bility of the field orientation of NACO over a period of more than 
7 months. From the analysis of archival NACO reference astro- 
metric fields, G. Chauvin (private communication) observed that 
a shift of 0.30+0.03 deg in the true North orientation occured be- 
tween April and August 2008. We identified an intervention on 
the instrument in early July 2008 that explains the observed shift 
(although it does not affect our April-May data). Other actions 
took place on the instrument later in 2008, that possibly modified 
or removed part of this rotation shift. The second cause is that the 
true accuracy of the orbital trajectory of B in our model (Fig. |3]l 



is limited by the accumulation of eiTors with time. Because of 
the uncertainties on the A-B orbital parameters, it is practically 
impossible to relate our 2009 measurement unequivocally to the 
2008 epochs at a level of a few 10 //as. We therefore excluded 
the isolated 2009 epoch from the present astrometric analysis. 

For the two epochs in early May 2008, however it is more 
difficult to reach a solid conclusion on the origin of the deviation. 
We considered the following possible explanations: 

• Atmospheric or instrumental effect. We checked the dif- 
ferent observational parameters of all the observations, either 
environmental (seeing conditions, airmass, etc.) or instrumen- 
tal (flat or bias properties, etc.), but we could not find a signif- 
icant particularity of the data collected over these two nights. 
The observed deviation was present over two nights separated 
by only two days, and was not present before or after these two 
nights. Although a peculiar instrumental bias cannot be formally 
excluded, this points to an astrophysical origin. 

• Transit of a dark body (exoplanet or faint star). The ge- 
ometry of the (5 Vel A-B eclipsing system makes it a favorable 
configuration for the transit of dark bodies in front of the stars of 
the eclipsing pair. The existence of these systems was predicted 
by [Sc hneider (1994), and recently demonstrated by Armstrong] 
|et al. (|2012j|. A search for transit signatures in the photometry 
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1.5 1.0 0.5 0.0 -0.5 

arcseconds 

Fig. 3. Relative orbit of (5 Vel B around ^ Vel A. The trajectory com- 
puted using the orbital elements determined in Paper III is represented 
as a black ellipse, and the measurements corresponding to our NACO 
and VISIR astrometry are shown in the enlarged insert box. The axes 
are labeled in arcseconds, and the numerical values of the positions are 
listed in Table|3] 

measured with SMEI (Fig.|6]l did not give positive results down 
to about 2% of the total Aa+Ab+B flux. There is still a possi- 
bility that a very shallow transit of 1-2% of the total flux may 
be present 4 - 5% of the flux of one of the main components 
Aa or Ab), but it would be insufficient to explain the observed 
astrometric signal. So the transit of a dark component in front of 
one of the stars appears unlikely. 

• Activity of 6 Vel Aa and/or Ab in the Bry line. The NACO 
images were taken with a narrow-band Bry filter. It is therefore 
possible that episodic mass loss created a temporary envelope 




2 -2 -1 1.0 0.5 0.0 -0.5 -1.0 

E - (mas] E ^ (mas] 



Fig. 4. Left: Visual rendering of the Aab system close to the primary 
eclipse (see also |Merand et al.|2011^ . Right: Observed astrometric dis- 
placement of the center of light of the Aab system from NACO differ- 
ential astrometry (points) with the model trajectory of the center of light 
(thick curve). The orbital trajectories of the Aa and Ab components are 
represented using thin curves. The rms residual dispersion of the mea- 
surements compared to the model is llOyuas (excluding the 2008 May 
05, 2008 May 07 and 2009 Ian 07 epochs, that are also not shown on 
this plot, see Sect. |4.1 [ for details). 

and circumstellar emission around one of the two eclipsing stars. 
An emission of this kind could be faint in the broad-band pho- 
tometry of the system (SMEI data), but would be much more 
apparent in our narrow-band Bry filter. This would affect our 
astrometric measurements by displacing the CL of the eclipsing 
system towards the star experiencing the circumstellar emission. 
According to our model of the system (Paper III), the expected 
Bry flux ratio between 5 Vel A and B is 10.67. We measure a 
median value of f(A)/f(B) = 10.37 ± 0.13 over all our obser- 
vations (Table |2]i, which is slightly below the model, but could 
be explained by interstellar extinction. It should be noted that 
the flux ratio values observed on 2008-05-05 and 05-07 are both 
lower than the median value (10.27 to 10.30). Moreover, the ob- 
served dispersion of the measurement points 0.13) is much 
larger than the uncertainty of the measurements (» 0.01). This 
means that either Aa, Ab, or B experiences photometric variabil- 
ity in the Bry line. The SMEI photometry presented in Fig.[6]also 
shows a possible variability at a 2-3% level. If we consider only 
the case of a change in brightness of one of the eclipsing stars (a 
change in the flux of B would have no effect on the astrometry), 
a 1% change in total flux means a 2% change in the relative 
brightness of 6 Vel Aa and Ab. This change in flux would result 
in a geometrical shift of the CL of the system of approximately 
300 /ias when the stars are at maximum elongation (separation 
16 mas). This amplitude is similar to the observed drift of the 
CL on 2008 May 05 and 2008 May 07. However, the fact that 
we have only the relative photometry of A and B at high preci- 
sion means that we cannot exclude a change in the brightness of 
B that would invalidate this hypothesis. 

• Binarity of 6 Vel B. The presence of an undetected com- 
panion around B would cause a displacement of the photocenter 
of this component, and thus aff'ect the differential A-B astrome- 
try. The mass of the companion of B would have to be relatively 
high, and its orbital period short, but we cannot formally exclude 
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Fig. 5. Top: Apparent displacement of the center of light of the 
6 Vel Aab eclipsing pair relative to B as a function of the orbital phase 
(the secular displacement of B has been subtracted). The two discrepant 
measurements around = 0.7 are discussed in Sect. |4. 1| Bottom: Ra- 
tio of the photometric fluxes of ^Vel A and B in the NACO narrow- 
band 2.17/jm filter (Table[2j, normalized to the median measured ratio 
(10.372). For both plots, the solid curves are the predictions of the 
model presented in Paper III, and not the result of a fit to the data. 



of Aa and Ab introduces an additive noise on the measurement 
of the Aab CL displacement. Although the formal accuracy of 
each individual astrometric measurement is unaffected, the level 
of agreement with the CL displacement predicted by our model 
(that assumes non-variable stars) is reduced. 

4.2. Astrometric accuracy 

The dispersion of our NACO astrometric measurements around 
our model of the system (110/vas rms, excluding the 2008 May 
05 and 2008 May 07 epochs, see Fig. |4]l is larger than our in- 
dividual measurement uncertainty (» 50;uas per epoch). We 
therefore estimate that the true accuracy of a single measure- 
ment epoch is most probably between 50 and lOOyuas. This 
level of accuracy is consistent with the conclusions of Fritz et| 



al. ( 2010| l. This confirms that very high-precision astrometric 



observations are feasible from the ground with adaptive optics 
systems on large telescopes. The measurements we presented 
were obtained on a particularly favorable object, very bright and 
with only two unresolved objects in the field of view. However, 
the detected astrometric displacement is particularly subtle, and 
the fact that we could retrieve the orientation of the orbital plane 
of the eclipsing pair demonstrates the potential of this technique. 

One should also note that we used a neutral density filter and 
a narrow band filter to prevent the saturation of the CONICA 
detector. This means that the level of accuracy we demonstrated 
is directly applicable to significantly fainter targets than 5 Vel 
(probably by at least 10 magnitudes). For even fainter objects 
(niK > 14), particularly in the Galactic Center region, a detailed 
discussion is presented by |Fritz et al.| ( |2010| . 



5. Conclusion 






53200 53400 
HJD-2400000 



Fig. 6. Total flux of the S Vel A-B system as a function of the phase 
from SMEI observations, excluding the eclipses, normalized to its aver- 
age over the observation period. 



this possibility. As suggested in Paper II, spatially resolved spec- 
troscopy of 6 Vel B would be necessary to check this hypothesis, 
as well as to confirm the physical properties of this star. 

In summary, the most likely explanations for the inconsistent 
astrometric measurements of 2008 May 05 and 2008 May 07 are 
either an unidentified instrumental bias, or a transient activity 
of one of the fast rotating Aa or Ab stars in the Bry hydrogen 
line (e.g. the temporary formation or disappearance of a disk 
contributing ^ 1 - 2% of the flux in the Bry line). The sensitiv- 
ity of our differential astrometric measurements to the flux ratio 



We reported high precision difl^erential astrometric measure- 
ments of the triple system 6 Vel. They show that it is possible 
to measure the separation of close binaries of suitable brightness 
with a repeatable precision of 50 to 100 //as. The comparison of 
the measured astrometric displacement of the CL of the eclipsing 
component of 6 Vel with a well-constrained model of the system 
shows that our NACO measurements give a realistic astrometric 
signal, in terms of amplitude and orientation on the sky. How- 
ever, this qualitatively good agreement is affected by what we 
interpret as photometric variability of one of the components of 
the eclipsing binary (although other explanations are also possi- 
ble). 

High accuracy astrometric measurements have the potential 
to bring novel constraints in several fields of astrophysics, from 
the Galactic center (e.g.|Gillessen et al.|2009j , to faint compan- 
ion detection (e.g. [Tokovinin et al.||2012|l, or cl uster dynamics 
(e.g. Clarkson et al.|20l"2f HuBmann et al.|2012[ l. The future in- 
strumentation of extremely large telescopes (see e.g. [Trippe et| 
laT 2(TlO) will give access to much fainter objects, although the 
actual measurement accuracy will probably not be improved sig- 
nificantly beyond the 50 juas limit. The next generation interfero- 
metric instrumentation at near-infrared wavelengths, particularly 
the GRAVI TY (Gil lessen et"aLl|2010i .E isenh aueret al.||201 1 



Vincent et a l. 201 1 ) and ASTRA lWoiUez et al.|2012| [Stone et] 



al.|2012 1 instruments, are expected to overcome this limit, down 
to an accuracy of 10 to 20 fias over a field of view of a few arc- 
seconds. 
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Table 3. Observed and model positions of (5 Vel B relative to A for the epochs presented in Fig.p](see also 



Argyle et al. 



20021. 



MJD 


obs(Q;) 


moael(a') 




obs(o) 


model((5) 




Aff 




Ao 


13467 


7500 


0.11 ±0.2 


0.29308 




-2.26 ± 0.2 


-2 


12963 


-0 


17874 


-0 


12748 


15100 


3000 


0.22 ± 0.2 


0.40581 




-2.30 ± 0.2 


-2 


33885 


-0 


18842 





03910 


18104 


7000 


0.48 ± 0.2 


0.59401 




-2.86 ± 0.2 


-2 


61761 


-0 


11537 


-0 


24262 


19753 


3600 


0.69 ±0.1 


0.68520 




-2.63 ±0.1 


-2 


71535 





00502 





08438 


21553 


8500 


0.81 ±0.1 


0.77407 




-2.89 ±0.1 


-2 


77964 





03269 


-0 


10985 


24132 


7000 


0.95 ±0.1 


0.88042 




-2.88 ±0.1 


-2 


79750 





06595 


-0 


08092 


25568 


0800 


0.85 ±0.1 


0.92825 




-2.67 ±0.1 


-2 


77070 


-0 


07696 





10325 


27996 


4500 


1.02 ±0.1 


0.98936 




-2.74 ±0.1 


-2 


66705 





03325 


-0 


06803 


31481 


5500 


0.97 ±0.1 


1.02966 




-2.32 ±0.1 


-2 


39249 


-0 


06104 





07692 


34342 


1300 


0.99 ±0.1 


1.01591 




-1.99 ±0.1 


-2 


05765 


-0 


02202 





06694 


34418 


8800 


0.92 ±0.1 


1.01491 




-2.00 ±0.1 


-2 


04731 


-0 


09296 





05202 


34378 


0000 


0.96 ±0.1 


1.01545 




-1.99 ±0.1 


-2 


05283 


-0 


05359 





06313 


43828 


7000 


0.565 ± 0.01 


0.57324 




-0.277 ± 0.01 


-0 


27270 


-0 


00838 


-0 


00402 


43862 


7000 


0.567 ± 0.01 


0.57025 




-0.262 ± 0.01 


-0 


26477 


-0 


00289 





00261 


48074 


5000 


0.116 ±0.01 


0.10948 




0.716 ±0.01 





69105 





00643 





02462 


48476 


0000 


0.058 ± 0.01 


0.05763 




0.755 ± 0.01 





76601 





00045 


-0 


01124 


48805 


0000 


0.008 ± 0.01 


0.01475 




0.774 ± 0.01 





82238 


-0 


00664 


-0 


04842 


51464 


4300 


-0.2886 ± 0.005 


-0.30349 


1 


0200 ± 0.005 


1 


00667 





01487 





01328 


54581 


5100 


-0.4259 ± 0.001 


-0.42983 





4541 ± 0.001 





45143 





00391 





00262 


54557 


0224 


-0.4306 ± 0.001 


-0.43010 





4578 ± 0.001 





45796 


-0 


00053 


-0 


00021 


54560 


9976 


-0.4307 ± 0.001 


-0.43005 





4571 ±0.001 





45690 


-0 


00068 





00020 


54562 


0121 


-0.4309 ± 0.001 


-0.43004 





4570 ± 0.001 





45663 


-0 


00082 





00033 


54563 


0048 


-0.4308 ± 0.001 


-0.43003 





4569 ± 0.001 





45637 


-0 


00073 





00051 


54576 


9715 


-0.4303 ± 0.001 


-0.42988 





4528 ± 0.001 





45264 


-0 


00044 





00013 


54579 


0231 


-0.4301 ± 0.001 


-0.42986 





4520 ± 0.001 





45209 


-0 


00019 


-0 


00012 


54580 


9917 


-0.4300 ± 0.001 


-0.42984 





4512 ±0.001 





45157 


-0 


00013 


-0 


00042 


54591 


9748 


-0.4290 ± 0.001 


-0.42972 





4485 ± 0.001 





44863 





00073 


-0 


00010 


54593 


9732 


-0.4294 ± 0.001 


-0.42969 





4481 ±0.001 





44809 





00030 


-0 


00004 


54604 


0442 


-0.4296 ± 0.001 


-0.42958 





4455 ± 0.001 





44540 


-0 


00005 





00010 


54838 


1347 


-0.4246 ± 0.001 


-0.42611 





3813 ±0.001 





38198 





00155 


-0 


00070 



Notes. All angles are expressed in arcseconds. The model positions were computed using the orbital elements in Paper III. 



